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Abstract 

Propagation of arbitrary-amplitude ion-acoustic solitary (IASWs) as well as periodic waves 
(IAPWs) is investigated in a fully degenerate quantum electron-ion plasma consisting of 
isothermal- or adiabatic-ion species. It is shown that the system dimensionality and degrees of 
freedom play critical roles in matching criteria for propagation of such waves. Furthermore, it 
is revealed that for the case of adiabatic-ion unlike isothermal one, in some cases, there exists 
an upper fractional ion-temperature limit for the existence of IAPWs. It is also shown that, the 
variations of wave-amplitude with respect to the change in fractional ion-temperature is quite 
different for the cases of isothermal and adiabatic-ion plasmas. 
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I. INTRODUCTION 



Ion acoustic waves (IAWs) have undeniable impact on scientific knowledge of character- 
istics of ionized environments. 



studied for long time since 1961 



inear and nonlinear features of these identities has been 
l| using diverse techniques, of those Sagdeev Q-Q and per- 
turbation |7H12| methods are well-known Recent achievements in micro- and nano-electronic 
device fabrication has motivated new efforts to study the interesting nonlinear behavior of 
IAWs in dense plasma environments such as quantum dots, quantum wires, quantum wells, 
carbon nano-tubes, quantum diodes, ultra-cold plasmas, micro plasmas, biophotonic, etc. 



131 ] . For instance, L. C. Gardner [14| has applied quantum hydrodynamic (QHD) model to 



study the electron-hole dynamics in semiconductors. 

Dense or degenerate warm plasmas are also abundant in astrophysical environments such 
as the core of planets, white dwarfs and neutron stars [15j]. It has been shown that in a dense 
plasma, where the quantum mechanical rules are in action, new nonlinear effects arise 



16- 



181 ] due to the non-classical quantum tunneling effect 19]. One of important parameters of 



a quantum plasma is the electron number- density which is a central parameter in definition 
of degenerate plasmas. As the number of electrons increase in quantum plasmas, quantum 
effects dominate due to Pauli exclusion principle leading to higher degeneracy pressure which 
is much larger than classical electron and ion pressures. Therefore, it is naturally expected 
that the system dimension and degrees of freedom to play crucial role in quantum plasma 
environments. 

In recent few years there has been intense investigations on nonlinear behavior of planar 
and non-planar IAWs in quantum plasmas due to their experimental as well as theoretical 
importance. Recent investigations indicate that the propagation characteristics of solitary 
waves in bounded quantu m p lasma with non-planar geometry significantly differ from those 



in unbounded planar one 



20j-|22|. Characteristic features of 



of IAWs has been studied extensively in electron-ion 23|, [24j as well as electron-positron-ion 



inear and nonlinear propagation 



quantum plasmas [25|, |26[. It has been found that only rarefactive large amplitude solitary 
excitations exist in quantum electron-ion plasmas, while, for electron-positron-ion plasmas 
only compressive large amplitude solitary excitations are allowed for given fractional plasma 
parameters. 

On the other hand, it seems that less attention has been paid to the system dimensionality 
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and ion degrees of freedom effects, which are supposed to be of vital importance in 
quantum systems. In current research we aim at exploring these effects on propagation of 
large-amplitude IAWs in an un-magnetized electron-ion quantum plasma. This investigation 
can be of sufficient importance in laboratory plasma research. The organization of the article 
is as follows. The basic normalized plasma equations are introduced in section [III Nonlinear 
arbitrary-amplitude solutions are derived in section III II Numerical analysis and discussions 
is given in section IIVI and final remarks are presented in section [V] 



II. THE BASIC EQUATIONS 

Lets consider a dense plasma consisting of inertialless degenerate electrons and iner- 
tial positive hot ions. It is assumed that the plasma is fully degenerate so that the ion 
temperature is much less than the characteristic Fermi-temperatures of electrons. It is also 
noted that in a fully degenerate plasma the rate of electron-ion collisions are limited due 
to Fermi blocking process, hence, the plasma is considered almost collision-less. Therefore, 
it is reasonable to use an isothermal pressure model for ions. Using the conventional fluid 
quantum hydrodynamics (QHD) basic equations, we have for both electrons and ions 



drij _|_ dnjVj 
dt ' dx 

dn e _|_ dn e v e 
dt i dx 



0. 



dvj , v .§Vi e dip 

dt 1 dx rrii dx 



( ■yk B T i \ (JH_\1- 2 drn I 
^m,n<») J \n(°)J dx " r 



7-2 ^ h 2 q 
2mi 2 dx 



d 2 /dx 2 [^/n-] 



(1) 



irii I dt c dx J rrti dx 



1 dP e 

mine, dx 



+ 



h 2 d 

2m e irii dx 



d 2 /dx 2 [^\ 



^ = 4vre (n e 



11; 



where 7 = 1 (7 ^ 1) with 7 = (/ + 2)/f (f being the ion degrees of freedom) corresponds 
to isothermal (adiabatic) ions case. Note that the left hand-side of momentum equation 
for electrons has been set to zero due to small electron-to-ion mass ratio. By using the 
following scalings 
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where, u P i = y 4ire 2 n^ /m>i and c s = \j2Ep e /mi are characteristic plasma-frequency and 
quantum sound-speed (Ep e = ksTpe is the electron Fermi-energy), respectively, we obtain 
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the normalized basic QHD equation set as 
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where the fractional temperature a (cr = measures the ion temperature rela- 

tive to the electron Fermi-temperature and is much less than unity in a fully-degenerate 
zero-temperature Fermi-gas. The new parameter, H = fihJ ve j1Ep e , is called the quantum 
diffraction parameter which is the ratio of plasmon-energy to electron Fermi-energy. Note 
also that, we have ignored the quantum term (last-term) in momentum equation of inertial 
heavy ions. 

One notices that, under the zero-temperature assumption in a dense and degenerate 
plasma environment, the degeneracy pressure follows from Pauli exclusion principle and 
relates to the particle number- density through the following relation 



( d+2 \ 



-Tip 



(4) 



e (d + 2) 

where, d represents the dimensionality of the system. Introducing the new variable £ = 
x — Mt (M is a measure of waves group-speed relative to that of quantum sound, c s ), and 
solving for the ion density and momentum equations altogether, we get 
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- oTnrij, (7 = 1), 
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Moreover, using Eq. ([41) , in the electrons energy relation, after integrating once with 
appropriate boundary conditions, we obtain 

H 2 d 2 
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It is evident that the conditions lim n* e 



1 and lim 1 



are met. 
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III. ARBITRARY-AMPLITUDE NONLINEAR WAVES 



In this section we will employ Sagdeev approach to find a suitable pseudo-potential which 
appropriately describes the propagation of ion-acoustic (periodic and solitary) waves and 
obtain some criteria for existence of such waves. This is done by evaluation of Eqs. (J5J) and 
(J6]) at local near-equilibrium state n e ~ rii = n. Using the above mentioned equations and 
defining new variable N = i/n yields 
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or equivalently, 
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The well-known energy integral reads as 
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It is clearly noticed that the potential and its first derivative vanish at n — 1. The possibility 
of ion-acoustic solitary waves require the following conditions to simultaneously met 
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n=l 



It is further required that for at least one either maximum or minimum nonzero n-value, 
we have U(n m ) = 0, so that for every value of n in the range n m > n > 1 or n m < n < 1, 
U (n) is negative. In such a condition we can obtain a potential minimum which describes 
the possibility of an IASW propagation. 

On the other hand, for existing a periodic IAW we require that 



U(n)\ n=1 
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and that, for at least one either maximum or minimum nonzero n-value, we have U(n m ) = 0, 
so that for every value of n in the range n m > n > 1 or n m < n < 1, U (n) is positive. In this 
situation we obtain a potential maximum describing the possibility of an IAPW propagation. 
Inspection of the second- derivative of pseudo-potential shows that 

d 2 U(n)/dn 2 \ n=1 < if M 2 -a<l/d 
d 2 U(n)/dn 2 \ n=1 > if M 2 -a>l/d 

(12) 

d 2 U{n)/dn 2 \ n=1 < if M 2 + 1 a<l/d 
d 2 U(n)/dn 2 | n= i > if M 2 + ja > 1/d 

It can be easily confirmed that for 7 = 1 we always have lirn U(n) = 2M 2 /H 2 > and 

n— >0 

lim U(n) = —00 independent of other plasma parameters. On the other hand, for 7 7^ 1 

n— >+oo 

we have lim U(n) = 2M 2 /H 2 > 0, independent of other plasma parameters and for f > d 

n-»0 

we obtain lim U(n) = —00, but for / = d we have lim U(n) = —00 if tr < l/(d + 2). 

n— >+oo n— >+oo 

It is therefore concluded that for IASW to exist we must have M 2 — a < 1/d for 7 = 1 and 
M 2 +ja<l/d for 7 / 1. On the other hand for the possibility of IAPWs for 7 = 1 we 
must have M 2 — a > 1/d and in the case of 7 7^ 1 with / > d we must have M 2 + 70" > 1 jd 
and for / = d we must have M 2 + 70" > 1/d and also a < l/(d + 2). 

The stationary wave solutions corresponding to this pseudo-potential which satisfis the 
mentioned boundary-conditions, read as 
f n dn 

{ - &=± y, T^m- (13) 

IV. NUMERICAL ANALYSIS 

Numerical evaluation of pseudo-potential profiles indicates that only rarefactive (n m < 1) 
large-amplitude solitary waves exist in this quantum electron-ion plasma for all system di- 
mensionality (d) and the ion degrees of freedom (/). This is in agreement with the previously 
reported results in Ref. j^] for quantum electron-ion plasma with the special case of 7 = 3 
in one dimension. In current numerical scheme we only plotted the figures for the special 
case of (point-like ions) d = f. In Figs. 1(a) (7 = 1) and 2(a) (7 7^ 1) the regions in 
M-a plane is shown, where IASWs or IAPWs can occur. The curves in these figures are 
M 2 — a = 1/d and M 2 + 70" = 1/d borders for 7 = 1 and 7 7^ 1, respectively, and the 
horizontal lines for 7 7^ 1 define the upper a-limits (a = l/(d + 2)) of IAPW existence for 



plasma with adiabatic point-like (d = f) ions. The left/right side of curves in Fig. 1(a) is 
the region where IASWs/IAPWs can exist. On the other hand, for 7 / 1 IAPWs stabil- 
ity region (M 2 + 70" > 1/d), shown in Fig. 2(a), is also bounded by an upper sigma limit 
(a < l/(d + 2)) for point-like ions d = f. Note that the a = (cold-ion) limits which 
match exactly for the cases of plasma with isothermal and adiabatic ions gives a maxi- 
mum/minimum (M m = l/y/d) for adiabatic/isothermal-ion plasma case. It is also clearly 
observed that for d — 1 in isothermal-ion case, contrary to that of adiabatic-ion, there is 
always a region (for all possible values of fractional ion-temperature, a) of supersonic soli- 
tary excitations, while for the shown range of a only subsonic IASWs exist for d = 2 and 
d = 3. However, for the case of adiabatic-ion plasma all solitary excitations lie in the sub- 
sonic region, for all dimensionality values. It is further remarked that the increase in the 
value of fractional ion-temperature increases/decreases the possible Mach- number range for 
all values of dimensionality for isothermal/adiabatic-ion case. It is observed that, as the 
dimension of the system is increased the possible Mach-number range for IASWs decreases 
for all values of o in both isothermal-ion and adiabatic-ion cases. 

Figures l(b-d) and 2(b-d) show the variations of the Sagdeev pseudo-potential profiles 
with respect to one-by-one variations of different fractional plasma parameters. For in- 
stance, it is concluded that for both 7 = 1 and 7 ^ 1 cases the shape of pseudo-potential 
is significantly affected by the change in the system dimensionality, quantum diffraction 
parameter and fractional ion-temperature for fixed other plasma parameters. Figures 3(a-c) 
and 4(a-c) presents the solitary wave profiles corresponding to the pseudo-potential profiles 
shown in Figs, l(b-d) and 2(b-d), respectively. From Figs. 3(a) and 4(a) it is clearly re- 
marked that for both 7 = 1 and 7 7^ 1 cases the IASW-amplitude decreases as the system 
dimension is increased, where, other plasma parameters are fixed. On the other hand, it 
is remarked from Figs. 3(b) and 4(b) that the soliton amplitude is unaffected by increase 
in the value of quantum diffraction parameter H, for all 7-values, however, its width is in- 
creased. This is a commonly reported feature of solitary excitations for many dense plasma 



environments 
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- rsj]. Moreover, it is concluded from Figs. 3(c) and 4(c) that increase 
in the value of fractional ion-temperature increases/decreases the ion-acoustic soliton am- 
plitude for isothermal/adiabatic-ion case. Figure 3(d) depicts the solitary profiles for sub- 
and supersonic cases (with d = 1 and 7 = 1) which indicates that the supersonic soliton 
is shorter and wider compared to the subsonic counterpart. Also, Fig. 4(d) depicts the 
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solitary profiles for isothermal (thick curve) and adiabatic (thin curve) cases (with similar 
other parameters) which indicates that the isothermal-ion solitary wave is slightly higher 
and narrower compared to that of adiabatic-ion. 

Figures 5(a) and 6(a) reveal that the increase in the Mach-number in the case of periodic 
ion-acoustic waves causes the wave amplitude to increase. A comparison has been made 
between pseudo-potential profiles in Figs. l(b-d)/2(b-d) for IASWs and those in Figs. 5(b- 
d)/6(b-d) for IAPWs, respectively. In addition to the large variation in pseudo-potential 
shape due to change in the system dimensionality for the case of IAPWs (compared to 
those of IASWs), it is observed that contrary to the case of IASWs the periodic ion-acoustic 
wave-amplitudes increase as the dimensionality increases. It is observed that the amplitude 
of IAPW and IASW are affected in a quite similar manner by the variations in the value 
of quantum diffraction parameter, H . A striking difference between IAPW and IASW for 
both 7 = 1 and 7 ^ 1 cases is the behavior of amplitude variation due to change in the 
value of fractional ion-temperature, which can be noticed by comparing Figs. l(d)/2(d) 
with 5(d) /6(d), respectively. It is evident that increase in the fractional ion-temperature 
increases/decreases the amplitude of IASWs for isothermal/adiabatic-ion case, while it de- 
creases/increases the amplitude of IAPWs for isothermal/adiabatic-ion case. 



V. CONCLUSION 



Using the Sagdeev potential approach, we investigated the large-amplitude ion-acoustic 
solitary as well as periodic wave propagations in a degenerate quantum plasma with isother- 
mal or adiabatic ion species in different plasma system dimensionality, taking into account 
the effect of ion-temperature and have shown that the system dimensionality and ion degrees 
of freedom has crucial impact on the propagation of such waves. In the studied range of frac- 
tional ion-temperature the one dimensional system with isothermal ions contains supersonic 
solitary propagations, while 2- and 3-D systems do not. For plasma with adiabatic ions only 
subsonic IASWs could be found for all dimensionality and for the chosen range of fractional 
ion-temperature. For the case of adiabatic-ion unlike isothermal one there is an additional 
upper fractional ion-temperature limit for the existence of IAPWs. It was also shown that, 
the variations of wave-amplitude with respect to the change in fractional ion-temperature is 
quite opposite for the cases of IASWs and IAPWs and for isothermal-ion and adiabatic-ion 
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cases. 
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FIGURE CAPTIONS 



Figure- 1 

(Color online) (a) The stability regions of arbitrary-amplitude IASWs and IAPWs is 
shown in M-a plane with degenerate electrons and isothermal ions. The region to the 
left /right of curves indicate the stability region for IASW/IAPW. The thickness of the 
curve increases as the dimension becomes higher, (b), (c), (d) Variations in Sagdeev pseudo- 
potential profiles of IASW with respect to different plasma parameters. The change in the 
dash-size appropriately represents the change in the varied parameters in each plot. 

Figure-2 

(Color online) (a) The stability regions of arbitrary-amplitude IASWs and IAPWs is 
shown in M-a plane with degenerate electrons and adiabatic ions. The region to the be- 
low/above of curves indicate the stability region for IASW/IAPW. The horizontal lines 
correspond to the upper a- limits of IAWs for the case of d = f (point-like ions). The 
thickness of the curve increases as the dimension becomes higher, (b), (c), (d) Variations 
in Sagdeev pseudo-potential profiles of IASW with respect to different plasma parameters. 
The change in the dash-size appropriately represents the change in the varied parameters in 
each plot. 

Figure-3 

(Color online) The variation of the solitary wave profiles for isothermal-ion with respect to 
plasma parameters represented in Figs, l(a-c). (d) Wave profiles for supersonic (thick) and 
subsonic (thin) solitary waves. The change in the curve-thickness appropriately represents 
the change in the varied parameters in each plot. 

Figure-4 

(Color online) The variation of the solitary wave profiles for adiabatic-ion case with 
respect to plasma parameters represented in Figs, l(a-c). (d) Wave profiles for isothermal 
(thick) and adiabatic (thin) solitary waves with similar other plasma parameters. The change 
in the curve-thickness appropriately represents the change in the varied parameters in each 
plot. 
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Figure-5 

(Color online) Variations in Sagdeev pseudo-potential profiles of IAPW for isothermal-ion 
case with respect to different plasma parameters. The change in the dash-size appropriately 
represents the change in the varied parameters in each plot. 

Figure-6 

(Color online) Variations in Sagdeev pseudo-potential profiles of IAPW for adiabatic-ion 
case with respect to different plasma parameters. The change in the dash-size appropriately 
represents the change in the varied parameters in each plot. 
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FIG. 1: 
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o-=0.1, H=0.1, M=0.5, d=l,2,3 (Isothermal) 
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cr=0.1, H=0.1, M=0.3, d=l,2,3 (Adiabatic) 



cr=0.1, H=0. 1,0.2,0.3, M=0.3, d=3 (Adiabatic) 



1.0 K 



0.8 



,0.6 



0.4 



0.2 




1.0 H 



0.9 



0.8 



0.7 



0.6 



0.5 



1.0 




cr=0,0.05,0.1, H=0.2, M=0.3, d=3 (Adiabatic) 



cr=0.1, H=0.1, d=l (Isothermal & Adiabatic) 




1.0 



0.8 



0.6 



0.4 



0.2 




FIG. 4: 



17 



o-=0.1, H=0.1, M=l. 1,1.2,1.3, d=3 (Isothermal) 
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